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Abstract — Different ways of modelling quantum control sys- 
tems, formulating control problems and solving the resulting 
problems are considered and compared. In particular, we 
compare the performance of geometric and optimal control, as 
well as iterative techniques for optimal control design versus 
local gradient optimization using a Lyapunov-type potential 
function for two problems of general interest: global control 
of qubits and entanglement generation in the form of Bell state 
preparation. 

I. Introduction 

Motivated in part by the rapid growth of nanofabrication 
techniques and nanotechnology, as well as a surge of interest 
in novel applications of quantum effects such as quantum 
information processing, the control of quantum systems is 
a subject that has received considerable attention recently. 
Many control strategies for quantum systems have been 
proposed, from open-loop strategies for Hamiltonian engi- 
neering involving usually off-line design of control fields 
based on a model of the system, assuming knowledge of the 
initial state and the control objective, to closed-loop strate- 
gies mostly based on conditional quantum trajectories and 
continuous weak measurements [1]. Some variants of open- 
loop Hamiltonian engineering have enjoyed considerable 
success in various experimental settings [2], [3] in the form 
of adaptive open-loop techniques such as direct laboratory 
optimization [4], [5]. However, despite the obvious practical 
importance of selecting the best strategy, very little work 
has been done in comparing different control strategies 
in terms of their effectiveness, efficiency and robustness. 
This is the topic of this paper. We will restrict ourselves 
here to comparing open-loop control strategies, and focus 
specifically on geometric, Lyapunov and iterative optimal 
control design techniques. 

There are many ways of formulating problems in quantum 
control as optimal control problems, from time-optimal con- 
trol [6] to variational techniques [7], [8], [9], [10] incorpo- 
rating various constraints, e.g., on pulse energies, amplitudes 
and frequencies. In addition, many quantum control problems 
can be formulated either as state control problems (using 
either a pure state vector, wave-function, density operator 
or real coherence vector representation), as optimization 
of observables problems, or as process control problems 
(e.g. producing a desired quantum logic gate). Unfortunately, 
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although these mathematical models are often equally valid, 
if not equivalent, ways of describing the system and control 
objective, we find that different formulations of the problem 
can lead to different solutions for the optimal control field. 
This is significant, as the solutions are not always equally 
desirable from a practical or physical point of view. Different 
formulations have different complexities when it comes to 
solving the resulting equations, either numerically or analyt- 
ically. 

Finally, once a particular model and problem formula- 
tion have been chosen, there are various ways of solving 
the resulting control problem, ranging from direct single- 
step control design methods based on Lyapunov potential 
functions [11], [12], to iterative techniques for control field 
design using gradient-based optimization algorithms such 
as GRAPE [10] and other monotonically convergent algo- 
rithms [7], [8], [9]. Unfortunately, again, even if the formu- 
lation of the control problem is the same, different techniques 
can yield very different solutions. The choice of algorithm for 
a particular control problem may depend on further factors. 
Certain algorithms are more flexible, allowing more realistic 
models to be used, such as the inclusion of non-Hamiltonian 
processes or a realistic restriction on the capabilities of the 
pulse-shaping equipment. Other algorithms are numerically 
and computationally far simpler — finding a global optimal 
solution is simpler for geometric control (where it is equiv- 
alent to finding a geodesic, and has an analytic solution for 
many systems) than for an optimal control problem (where 
the numerical solution and discretization of the solution is 
non-trivial to globally optimize). The techniques also have 
different degrees of robustness to errors in our models of 
the systems. These competing issues of numerical accuracy, 
facility of implementation and robustness to modelling and 
systematic error must be appreciated to find the best practical 
solution. 

Since we obviously cannot address all of these problems 
in a single paper, we will illustrate some of the issues 
considering two typical problems of significance in quan- 
tum computing scenarios: simultaneous, selective control of 
several qubits using global control pulses, and controlled 
entanglement of coupled qubits. We will compare geometric 
and optimal control design, based on direct optimization of a 
Lyapunov function and iterative optimal control design using 
monotonically convergent algorithms. Before we delve into 
comparing different control design techniques in Sections |IV] 
and [V] we briefly consider different ways of modelling 
quantum control systems and formulating quantum control 
problems in Sec III] and review the relevant control design 
techniques in Sec. |lll| 



II. Modelling quantum control systems & 

CONTROL PROBLEMS 

As eluded to in the introduction, quantum systems can 
be modelled in various ways. Systems initially in a pure 
state, i.e., unentangled with their environment and governed 
by Hamiltonian evolution can be described by a wavefunc- 
tion encoding information about the state of the system 
that evolves according to a control-dependent Schrodinger 
equation 



zn|*(x,t)) = i/[x,f(i)]|vl/(x,t)) 



(1) 



where H[x.,i{t)] is the control-dependent Hamiltonian and 
f the control field applied. For convenience we will choose 
units here so that /i ~ 1. In general, the wavefunction of 
the system is a normalizable and differentiable complex- 
valued function of time t and space x, and H is a partial 
differential operator, but Eq. ([TJ can be recast as an ordinary 
differential equation, e.g., by expanding |^'(x, t)) in terms of 
the eigenstates \n) of Hq, the system's intrinsic Hamiltonian, 
|\['(x,t)) = ^^71=1 '^ri\n) , where N is the dimension of the 
Hilbert space TC, in which case we simply obtain a linear 
matrix ODE 



c{t) = -iH[i{t)]c{t), 



(2) 



where c is a complex column vector, usually normalized to 
unity, ||c|| = 1, and H[i{t)] is simply a control-dependent 
NxN Hermitian matrix. The dependence of the Hamiltonian 
on the control fields is furthermore often assumed to be linear 
so that 



M 



H[{{t)]= Ho+Y. frr.{t)Hrn, M<00. 



(3) 



This the standard model for pure-state quantum control. 

We can also represent the state of the system by an 
operator p{t) acting on the Hilbert space H. For a pure- 
state system represented by |^'(x, t)) = c{t), we have simply 
p{t) = |\I'(x, i)) (^'(x, t) I = cc''', where c''' is the Hermitian 
conjugate of c, from which we can immediately deduce the 
evolution equation 



(4) 



where [A, B] — AB — BA is the matrix commutator. For a 
pure-state system subject to Hamiltonian evolution the latter 
formulation is equivalent to Eq. (|2]). Its main advantage is that 
it can be extended to systems initially entangled with their 
environment, i.e., in a non-pure quantum state, or can interact 
(non-coherently) with their environment by adding a non- 
Hermitian super-operator L£)[/5(<)] to the RHS of Eq. 

Finally, we can represent the state of a quantum system 
in terms of the (real) expectation values of a complete set of 
(Hermitian) observables. There are obviously many choices 
for the set of observables but a very convenient choice are 



the normaUzed Pauli matrices 
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(-|r)(.| + |.)(r|) 

\fc=i 



r\r + l){r + l\ 



for 1 < r < — 1 and r < s < A^. It is easy to show that 
the iV^ — 1 trace-zero matrices ak above, together with the 
normalized identity operator (Tq = ^7^^^ form a complete 
orthonormal basis for the Hermitian operators on Ti,, and 

hence we have p = I]«=o ^ = ('^'fe I P) = 

Tv{cr\p). Thus, we can represent any quantum state p by a 
real iV^ — 1 vector s = (sfc)^|j~^, where sq = 1/VN is 
omitted as it is constant. The vector s is generally called 
the Bloch or coherence vector in the physical literature, 
or the Stokes tensor in the mathematical literature. This 
Bloch representation is very popular especially for two-level 
systems as it permits visualization of quantum states aspoints 
inside a ball in M^. [^It is very easy to see that Eq. ^ thus 
gives rise to the so-called Bloch equation 

s(<) = A[i{tMt), (5) 

and Eq. ^ implies A ^ Aq + Y.m=i /(*)^™' where A™ 
are real-orthogonal matrices given by Am — ad^ for 
m — 0, . . . , M, for Hamiltonian systems. For dissipative 
systems one can derive an affine-linear Bloch equation 
s{t) = A[i{t)]s{t) + b, although we shall not discuss this 
case here. 

As this brief overview shows, there are several ways of 
modelling quantum systems, which are often equally valid 
or even mathematically equivalent descriptions of the system. 
We shall see, however, that they are not necessarily equally 
desirable from a practical point of view, and depending on the 
problem considered, some may in fact be highly preferable 
to others. 

But even if we have chosen a particular way to model the 
system, most quantum control problems can be formulated 
in many different ways. For instance, the three 'canonical' 
control problems 

• State control: Steer the system from a given initial state 
to a desired target state 

• Process control: Implement a desired quantum process 
(unitary operator for Hamiltonian systems) 

• Observable control: Drive the system such as to opti- 
mize the expectation value of an observable 

are inter-convertible. E.g., a pure state control problem can 
be formulated as observable optimization problem where 
the observable is the projector onto the desired pure state. 
A mixed-state control problem can be formulated as an 
observable optimization problem with A = p. An observable 

'in the case Af = 2, it is customary to normalize the Pauli matrices such 
that (Tj; CTf = 2(5fef in order to obtain the unit ball in M"^ , but this is a minor 
detail. 



optimization problem can be translated into the problem of 
simultaneously driving a set of basis states to coincide with 
the desired eigenstates of the observable, etc. In addition to 
many obvious choices, there are often less obvious ones that 
may turn out to be easier to solve than the more obvious 
choices. 

III. Control design techniques for quantum 

SYSTEMS 

Before we compare control design techniques, we briefly 
review the relevant techniques. The simplest control design 
strategy, with a long history of successful application in nu- 
clear magnetic resonance experiments, is the use of geomet- 
ric control pulse sequences. Although it is rarely formulated 
this way in NMR spectroscopy, the basic idea of geometric 
control is to find a target (unitary) operator that will achieve 
the control objective — be that the implementation of a par- 
ticular process, or the preparation of a quantum state, etc — 
and to decompose or factorize it, using standard techniques 
such as the generalized Euler angle or Cartan decomposition, 
into elementary operations that can be practically realized by 
applying simple geometric control pulses such as piecewise 
constant controls or Gaussian pulses, for example. Typically, 
geometric control sequences can be visualized as performing 
a sequence of rotations on Bloch vectors in whence 
the name geometric control. 

An alternative to simple geometric pulse sequences is the 
use of temporally or spectrally shaped pulses optimized for 
a particular task. The design of optimally shaped pulses is 
a non-trivial problem, which is usually solved numerically 
using one of a number of (similar) iterative techniques [7], 
[8], [9], [10]. The particular algorithm we will employ in the 
following involves iteratively solving an initial value problem 
for a variational trial function Sv^\t) representing the state 
of the system 



= L[f(")(t), 4")(t)], 5(")(io) = So, (6) 

followed by a final value problem for a variational trial 
function Ay{t) (usually representing an observable) 

= L[r(t),4")(i)], A^:^\tF) = Af, (7) 

while updating the control field in each step according to 

A O/m 

starting with an initial trial field f'^''^(t), and setting 



Sf„ 



{Ay{t)\df^^tMi{t),Sym. (8) 



super-operator L[i{t), Sy{t)] depends on the chosen model. 
For a Hamiltonian pure-state system 

L[f{t),Sy{t)] = ~iH[i{tmy{t)), 

for a general mixed-state system 



L[f(i),^,(i)] 



H[{{t)],pyit) +LD[Pvit)], 



and for a unitary operator control problem --iH[i{t)]Uy{t), 
etc. For a system with control-linear Hamiltonian ^ the 
derivatives of the dynamic operator with respect to f„i{t) 
are explicitly 

I 5/„L[f(t),5,(t)]) = (x.W|i?m|*.(i)) (9) 

for a pure-state system, 

(A„(t) I dfMHt),Svm = Tr (Ay{t)[H,n,Pv{t)] 



In the abstract notation above, Sy could be a pure state 
I'^v), a density matrix py, a real Bloch vector s„, or even a 
unitary operator Ay is a suitable conjugate variable. The 



(10) 

for a mixed-state system, etc. a, (3, A and the target time 
tp are non-negative real parameters with < a,/3 < 2 to 
ensure convergence. It is possible to let a, (3 and A be positive 
functions rather than constants but we shall assume constant 
values in this paper 

The third technique we consider is model-based (open- 
loop) control design using Lyapunov functions. This ap- 
proach essentially involves choosing the control field at every 
point in time such as to ensure the dynamical evolution 
results in the monotonic decrease of a Lyapunov potential 
V, which has a minimum at our target state. If the objective 
is to reach a given target state, an obvious choice for y is a 
monotonic function of the distance between the current and 
target states of the system, e.g., V{p,pd) — \\p~ pdW^, if we 
use the density operator formulation, V^d) = llV'^V'dlP 
in the pure state formulation, or y(s,S(i) = ||s — s^jp 
if we use the Bloch vector representation of the state. 
Regardless of the precise formulation, one easily obtains a 
simple rule for choosing the control field. For instance, in 
the Bloch representation, given a system satisfying s{t) = 
A[J{t)]s{t) with A[f{t)] = ^0 + f{t)Ai where Aq and 
Ai are dynamical generators, we can easily obtain the rule 
f{t) = KSrf(i)^Ais(i), where k is a positive constant. The 
obvious advantage of this approach is that we obtain an 
explicit control law, and under certain conditions, asymptotic 
convergence from a near-global set of target states to the 
target state can be proved [11], [12]. 

IV. Geometric vs optimal control 

As an explicit problem in this section we consider an 
ensemble of five quantum dots, which will be modelled as 
a two-level systems (qubits) with slightly different energy 
levels. Obviously, this is a simplification. Real systems will 
have more energy levels but it is a good model system 
to study the limitations of frequency selective geometric 
control and explore what improvements might be gained 
using optimally shaped pulses. The control objective is 
simultaneous, selective control of all qubits using global 
control pulses, i.e., the ability to simultaneously perform 




Fig. 1. (a) Two Gaussian wave-packets with 'pulse area tt' with frequencies 
oJi (pulse 1, black envelope) and 0)3 (pulse 2, red envelope), respectively, 
(b) Optimally shaped pulse to achieve simultaneous selective excitation of 
dots 1 and 3 while minimizing the (final) excitation of all other dots. 

selective gate operations on several qubits without any local 
addressing. For the first example we consider local operations 
and assume negligible inter-qubit coupling on the time-scales 
relevant for single qubit gates. 

If the individual qubits have different resonance frequen- 
cies, due to variations in size, shape or composition of the 
dots (as fabricated systems such as quantum dots usually 
do), the simplest control approach is to employ frequency- 
selective geometric control pulses to selectively induce de- 
sired single-qubit rotations without local addressing. This 
approach is attractive due to its simplicity, and frequency- 
selective geometric control pulse sequences have been used 
successfully in nuclear magnetic resonance (NMR) appli- 
cations, and to a lesser extent in electron spin resonance 
experiments. However, in practice such simple control pulse 
sequences are generally not optimal (for gate operation times 
and other figures of merit), and problems are expected to 
arise, e.g., when the frequency separation between the indi- 
vidual qubits is small, and fast control operation is desired 
to minimize gate operation times and deleterious effects of 
decoherence. In this case, we expect optimal control designs 
to out-perform the simple geometric schemes. Indeed, sim- 
ulations show that this is generally the case. 

To consider a concrete example, we compare the Bloch 
trajectories for a simple geometric pulse sequence and an 
optimally shaped pulse for the problem of performing si- 
multaneous bit flips on two qubits (here 1 and 3) without 
affecting the others. (The bit flip operation was chosen as it is 
easy to visualize on the Bloch sphere. Similar results hold for 
other single qubit gates.) The geometric control solution in 
this case is simple: apply a sequence of two tt pulses resonant 
with the transition frequencies of the two target dots. If the 
control pulses are sufficiently long, the geometric control 
pulse sequence produces very good results. However, as we 
reduce the pulse lengths and increase the pulse strengths, 
off-resonant excitation becomes a problem, as Fig. [T] shows. 
The left column (a) shows the ideal evolution of each dot in 
the rotating frame subject to the fields in Fig [T] (a). Target 
dots 1 and 3 follow a smooth path from south to north pole, 
the others are unaffected. However, the middle column (b), 
which shows the actual evolution of the quantum dots in the 
stationary lab frame when subjected to the control fields in 
Fig [T] (a), shows clearly that the applied pulses induce far 



Ideal Trajectories Real Trajectories Real Trajectories 
(Rotating Frame) (Lab Frame) Optimal Field 




Fig. 2. Bloch sphere trajectories for five quantum dots (with resonance 
frequencies of 1.32, 1.35, 1.375, 1.38, 1.397 eV, respectively) subject to 
two Gaussian pi pulses of approximately 2 pico-seconds each (a,b), and 
a shaped pulse designed using iterative optimal control (c). The actual 
trajectories for the geometric control pulse sequence (b) differ significantly 
from the ideal model trajectories (a). 

more complex evolution. Dots 1 and 3 are resonantly excited, 
leading to population transfer from the south to the north 
pole along a spiral path in the stationary frame, as desired. 
However, there is also significant off-resonant excitation of 
the remaining dots, all of which are left in various excited 
states at the final time. The right column (c) shows the path of 
the dots subject to field [T](b): despite following complicated 
trajectories both target dots finish at the north pole and all 
the other dots are returned to the south pole (ground state) 
at the final time, leaving no unwanted excitations. 

The optimal pulse shown in Fig [T] (b) was designed by 
formulating the problem as an optimization problem for the 
observable A = diag(0, 1, 0, -1, 0, 1, 0, -1, 0, -1), defined 
on the Hilbert space H = ®\^iTik, where Ti^. is the 
single qubit Hilbert space. Although for five qubits the 
problem of directly optimizing the average gate fidelity on 
the 32-dimensional direct product space (Xif.^j^Tifc is com- 
putationally still tractable on modern computers, our rather 
unintuitive formulation was considerably more efficient here. 
Since we assumed inter-dot coupling to be negligible for 
pico-second pulses in our example, it is sufficient to model 
the evolution of the individual dots (for the purpose of imple- 



meriting single qubit gates) on the direct sum Hilbert space 
H = ©l^iTifc, which has only dimension 10. Furthermore, 
the observable A here assumes its maximum exactly when 
dots 1 and 3 are in the excited state |1) and all others are 
in the ground state |0). Since all our dots are initially in 
the ground state |0), optimizing (A) requires implementing 
a unitary operator U such that U\0)k = for fc = 1,3 
and U\0)k = \0)k for fc = 2,4, 5. For this simple example it 
is easy to see that the only unitary operator (in ©|^i6il(2)) 
that achieves this (on the direct sum Hilbert space) is indeed 

U = XQ)iS)XQ)iQ)i where X = ^ is the bit flip 

operator — the operator we wanted to implement. Similar ob- 
servables can be constructed to implement simultaneous local 
operations on N uncoupled qubits, although the approach 
does not straightforwardly generalize to qudit systems. 

The optimal control problem thus defined was solved 
using the iterative optimal control design algorithm described 
in Sec. [Ill] The Liouville operator for the system under 
consideration was 



L[fit),p{t)]^[Ho + fit)H,,pit)h 
r(fc) 



(11) 



where H,n = diag(i7,Vi ), fc = 1, . . . , 5, to = 0, 1, and single 
qubit Hamiltonians 



Mk)_ 0\ Mk)_ 1 / 



1 - i 




The optimal pulse shown in Fig. [T] (b) was obtained specif- 
ically for f^°\t) = 0.559 sin(i) with a = /3 = 1 and 
A = 4. Convergence for this choice of control parameters 
was comparatively slow (500 iterations), but it was found 
that choices for a, /3, A that resulted in faster convergence 
tended to produce pulses with undesirable features such as 
high frequencies or high amplitudes, or lower gate fideli- 
ties. The relationship between a,/3, A and the initial trial 
field, and the solution obtained, in particular its suitability 
for implementation, versus algorithm convergence time are 
complex, important and worthy of further investigation. 

V. Optimal control: Lyapunov vs iterative 

CONTROL DESIGN 

In the previous section we showed that optimal control 
design can yield significant improvements over geometric 
control for certain problems, and that unconventional formu- 
lations of a control problem may significantly reduce the 
complexity of a problem and improve the quality of the 
results. In this section, we will compare two basic strategies 
for optimal control design: the iterative approach that was 
employed in the previous section and the mathematically 
more elegant and computationally less expensive approach 
based on Lyapunov functions. 

Although we could consider the same problem as in the 
previous section, for variety we shall instead consider the 
problem of preparing a Bell state l^f^) = ^(|00) + |11)) 
starting with a product state |00) for two weakly-coupled 
spins evolving under the control-dependent Hamiltonian 



H 



[fit)] - fit) (0.941) + 42)) + 0.1a, ® a„ (12) 



a typical (model) Hamiltonian for systems involving nuclear 
spins, e.g., in NMR. In this case we formulated the problem 
as a straightforward state control problem. Fig. [3] (top) 
shows the results produced by the iterative optimal control 
algorithm for A chosen to be the projection onto the target 
state and a = /3 = A = 1 for a target pulse duration of tp = 
200 time units. While the coiTect choice of the parameters 
a, (3, A, /*^°), ... in the iterative scheme appears to be crucial 
to obtain good results, and finding suitable parameters is 
not always a trivial task, in most cases it was obvious very 
quickly when we had chosen the parameters poorly, allowing 
us to adjust the parameters. After some tuning, the control 
pulses obtained from the iterative scheme were efficient in 
steering the system very close to the target state and had 
desirable characteristics, as the example shown in Fig. [3] 
shows, despite the fact that system (12i is not controllable. 
The distance from the target state at the specified target time 
is very small, and the plot of the populations and relevant 
coherence |pi4| shows smooth trajectories approaching the 
target values almost monotonically. Thus, the pulse is highly 
efficient in steering the system to the target state, near energy 
optimal, and has desirable characteristics such as limited 
variation in amplitude (no spikes) and a narrow frequency 
spectrum centered around the transition frequencies of the 
system. 

The Lyapunov potential approach proved considerably 
more problematic. Firstly, many initial states of interest 
such as the separable state |00) in our example, lie within 
the set of states for which convergence does not occur as 
(^'d(0)|iJi|^'(0)) = 0, and hence the Lyapunov derived 
control field / = 0. We can overcome this problem by 
applying a random initial pulse to kick the system out of 
its initial equilibrium state. Even with this modification, 
however, in all our experiments the controls derived from the 
Lyapunov potential function performed significantly worse 
than those derived from iterative schemes, both in terms 
of the rate of convergence to the target state, and the 
characteristics of the pulses. In fact, in our example, the 
distance from the target trajectory, although monotonically 
decreasing as a function of time, exhibits plateaux where 
the system appears to be trapped in 'false' semi-stable states 
for long periods of time. Furthermore, the control fields we 
obtained, a typical example of which is shown in Fig. [3] 
(bottom), tended to be generally spiky and spectrally more 
complex than those derived from iterative schemes. 

In our particular example (and a number of other cases 
we have studied so far) we found that the iterative algorithm 
allowed us to find control pulses that steer the system to 
a state very close to the target state in a relatively short 
time, as Fig. |3] shows. This applied even to some systems 
that are not generally controllable, provided that the control 
objective was dynamically realizable, and in cases where it 
was not attainable, we were often able to steer the system to 
reachable states that maximized the expectation value of the 
target observable subject to the dynamical constaints. 

The local optimization approach on the other hand, despite 
recent claims of strong global convergence [12], [13], seemed 



After 1 00 iterations (a=1 .00, |3=1 .00, X=^ .00) 
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Fig. 3. Control pulse, con'esponding evolution of the populations pi^i^, 
k = 1,2,3,4, and crucial coherence |pi4|, and distance from the target 
state as a function of time for an optimal control pulse obtained from an 
iterative scheme (top) and a Lyapunov-derived pulse (bottom). Insets shows 
the main Fourier components of the pulse. 

to fare significantly worse in our simulations. The system 
often got trapped for long periods in semi-stable critical 
points, and convergence to the target state was slow at best. 
Thus, while asymptotic convergence may be an attractive 
property mathematically, it appears to be less useful in prac- 
tice, especially where gate-operation and state preparation 
times are critical, although further work is necessary to 
assess whether the cases we have studied are anomalous, and 
whether the performance of the local gradient optimization 
can be improved, e.g., by choosing different Lyapunov poten- 
tials. The latter question is significant as local optimization 
techniques can in principle be adapted for measurement- 
based based feedback control in the context of stochastic 
differential equations [14], [15]. 

VI. Conclusion 

We have considered different ways of modelling quantum 
control systems, formulating open-loop control problems 
and techniques for solving them. The preUminary results 



presented here strongly suggest that model-based optimal 
control design tends to produce control pulses that are more 
effective, robust and efficient, than simpler schemes based on 
geometric pulses. However, the quality of the pulses obtained 
can vary significantly depending on the choice of model, 
problem formulation and the technique used for control 
design. 

We have seen in particular that control design based on 
instantaneous minimization of a Lyapunov potential function, 
despite attractive features such as simplicity and asymp- 
totic convergence properties, tends to produce fields inferior 
to those obtained from common iterative optimal control 
schemes in various regards, from rate of convergence to 
robustness and what might be termed energy and spectral 
optimality. Preliminary results suggest that iterative optimal 
control design produces far better results. 

Nonetheless, instantaneous optimization techniques re- 
main of interest due to the potential to adapt these techniques 
in the context of model-based feedback control. More work 
is necessary to understand how model choice, problem 
formulation and algorithmic parameters affect the conver- 
gence behavior and characteristics of the fields obtained, in 
particular for instantaneous optimization techniques, and how 
these characteristics might be improved, e.g., by choosing 
more sophisticated Lyapunov potentials. 
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